Thermalization of magnetically trapped metastable helium 
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We have observed thermalization by elastic collisions of magnetically trapped metastable helium 
atoms. Our method directly samples the reconstruction of a thermal energy distribution after the 
application of an RF knife. The relaxation time of our sample towards equilibrium gives an elastic 
collision rate constant of a ~ 5 x 10~^cm''/s at a temperature of 1 mK. This value is close to the 
unitarity limit. 



Bose-Einstein condensation (BEC) of dilute atomic va- 
pors has been observed in Rb g, Na §, Li § and H [|. 
Atoms in these gases are in their electronic ground state. 
Metastable helium in the 2 ^Si state (He*), which has 
long been of interest to the laser cooling community, is 
by contrast in a state 20 eV above the ground state. This 
situation presents new possibilities for the study of cold 
dilute atomic gases. First, the large internal energy per- 
mits efficient detection by ionization of other atoms and 
surfaces: it is possible to study very small samples. Sec- 
ond, Penning ionization by both the background gas and 
between trapped atoms offers a high time resolution mon- 
itor of the number and density of trapped atoms. Third, 
the possibility of using the large internal energy of He* 
for atomic lithography has already been demonstrated 
1^, and this application as well as atom holography ||] 
may benefit from highly coherent sources. Finally, much 
theoretical work has already been devoted to estimation 
of the elastic collision cross sections on the one hand and 
Penning ionizing rates on the other |7[ Experiments 
such as the one reported here can test this work. 

BEC is achieved in dilute gases by evaporative cool- 
ing of a magnetically trapped sample In He*, it is 
hampered by the fact that in a magneto-optical trap, the 
typical starting point of magnetic trapping, the achiev- 
able atomic density is limited by a large light-assisted 
Penning ionization rate |ll|, ^ |l^, |l5|. On the 
other hand, the scattering length for low energy elastic 
collisions is predicted to be quite large, and the Penning 
ionization rate highly suppressed in a spin polarized sam- 
ple H . He* in a magnetic trap necessarily constitutes 
a spin polarized sample and experiments have already 
demonstrated a suppression of more than one order of 
magnitude |8j. If the theoretical estimates are 

right, efficient evaporative cooling may still be possible in 
spite of the low initial trap density. We report here the 
observation of the thermalization of He* due to elastic 
collisions which appears to roughly bear out the predic- 
tions. 

To perform a thermalization experiment, a trapped 
cloud is deliberately placed out of equilibrium and its 
relaxation due to the elastic collisions between trapped 
particles is observed. Usually the observations are made 
by imaging the spatial distribution as a function of time 



|]l9| , po| |2l|| . In our experiment the relaxation is observed 
in the energy distribution of the atoms in the magnetic 
trap. First this distribution is truncated above -Erf = hv 
by a radio-frequency pulse (or RF knife) of frequency v. 
The cloud rethermalizes by elastic collisions and the pop- 
ulation of the states of energy higher than E^if increases 
from zero; for large times compared to the thermalization 
time Tth the distribution reaches a thermal distribution 
p2t . With the help of an analytical model and numerical 
simulations, we deduce Tth from the time dependence of 
the number of atoms with energy above E^rf- We mea- 
sure this time dependence by applying a second RF knife 
after a delay time and with a frequency slightly above 
that of the first one. Our model also allows us to relate 
Tth to the elastic collision rate per atom in the trap. 

Much of our setup has been described previously 
pTj . Briefly, we use a LN2 cooled DC discharge source 
to produce a beam of metastable He atoms. The beam 
is slowed down to 100 m/s using Zeeman slowing and 
loads a magneto-optical trap. Typically, 3 x 10^ atoms are 
trapped at a peak density of 3 x 10^ at/cm^^, limited by 
light-induced Penning ionization. The temperature of the 
cloud is about 1 mK and the cloud is roughly spherical 
with an RMS size of 2.5 mm. We then apply a 5 ms 
Doppler molasses to cool the atoms down to 300 //K. This 
is achieved by switching off the magnetic field, decreasing 
the detuning close to resonance and lowering the intensity 
to 10 % of its value in the MOT. An optical pumping 
step allows us to trap up to 1.5 x 10^ atoms in a loffe- 
Pritchard trap. We use a "cloverleaf configuration |2^ 
with B' = 85 G/cm, B" = 25 G/cm^ and a bias field 
Bq = 200 G. The two sets of coils are outside the vacuum, 
separated by 4 cm. After lowering the bias field to 4 G, 
the temperature of the compressed atomic sample reaches 
1 mK. The lifetime of the trap is 60 s. 

We use a 2 stage microchannel plate (MCP) to detect 
the atoms. The MCP is placed 5 cm below the trap- 
ping region and has an active area of 1.4 cm diameter. 
Two grids above the MCP allow us to repel all charged 
particles and detect only the He*. After turning off the 
magnetic trap, the MCP signal corresponds to a time of 
flight spectrum (TOF) which gives the temperature of 
the atoms. The area of this spectrum is proportional to 
the number of atoms in the trap at the time it was turned 
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FIG. 1: RF spectrum of atoms in the magnetic trap, a) Fraction of remaining trapped atoms after the RF pulse as a function 
of the RF frequency u. b) Derivative of these data, i.e. the energy distribution in the magnetic trap. The solid line is the 
prediction for a cloud at a temperature of 1.1 mK, the temperature measured by time of flight (TOF). The dashed line indicates 
the frequency corresponding to the bias field. 



off. The collection and detection efficiency of the MCP 
varies by roughly a factor of two depending on the mag- 
netic field configuration we use, and so one must take 
care to only use data corresponding to the same mag- 
netic field when making comparisons. We also use the 
MCP to monitor the atoms falling out of the trap while 
applying an RF knife. The area of the MCP signal in this 
case measures the number of atoms with an energy above 
that of the RF knife. Finally, when we bias the grids so 
as to attract positive ions, the MCP signal can be used 
to observe the products of Penning ionization with the 
background gas while the trap is on. This signal is pro- 
portional to the number of trapped atoms. We observe 
an exponential decay, indicating that two body loss (He* 
-I- He*) is neghgible. 

Two parallel coils in the vacuum system produce an RF 
magnetic field perpendicular to the bias field and consti- 
tutes the RF knife. To understand the effect of the RF 
knife on the trapped cloud and to assure that our sam- 
ple is at thermal equilibrium, we first performed an RF 
spectroscopy measurement of the energy of the atoms in 
the trap We apply an RF pulse at a frequency hv 
which changes the Zeeman sublevel of the atoms from the 
trapped M — +1 state to M = 0. The duration of the 
knife is 3 s, which is necessary to expel all the atoms with 
energy above hv over the entire range which we explore. 
We then turn off the magnetic trap to measure the num- 
ber of remaining atoms. Observation of the atoms falling 
onto the MCP during the RF knife shows that the flux of 
atoms expelled is negligible at the end of the pulse. An 
example of the RF spectrum is shown in Fig. |^a). The 
derivative of the data gives the energy distribution. In 
Fig. |l|b) we compare this distribution with a thermal one 
at 1.1 mK, the temperature measured by an independent 
TOF measurement. We conclude that our atomic sample 
is close to thermal equilibrium. 

We begin the thermalization experiment with a 2 s 



RF knife of frequency vi = 135 MHz (corresponding to 
■q — {hv — 2iJ,BBo)/kBT ~ 6). Next we measure the 
number of atoms falling onto the MCP during a second 
RF knife at a slightly higher frequency {1/2 = 138 MHz) 
and delayed by a time t. Assuming that the angular dis- 
tribution of the atoms expelled by the second RF knife 
is constant during the thermalization process, the MCP 
signal is proportional to the number of expelled atoms. 
Plots of the number of expelled atoms as a function of t 
are shown in Fig. ^ for samples having different numbers 
of atoms but the same temperature to within 10%. Fig. ^ 
shows that the number of atoms above the RF knife in- 
creases rapidly and then falls again with a time constant 
close to the trap lifetime as atoms are lost. If the ini- 
tial increase is indeed due to thermalizing collisions, the 
initial slope of each curve should be proportional to the 
square of the number of atoms. Our data roughly confirm 
this dependence. 

To be more quantitative, and to determine the ther- 
malization time Tth, we use a model based on the Boltz- 
mann equation under the sufficient ergodicity hypothesis 
and inspired by ||25|| . We divide the sample into two en- 
ergy regions, and with energies below and above 
ryfeeT respectively and denote by N- and A^-)_ the num- 
ber of atoms belonging to the two regions. We assume 
that 77 ^ 1. Immediately after truncation, = 0, and 
we seek a differential equation governing the time depen- 
dence of N+. Since 77 3> 1, we only take into account 
collisions of the type {£-) + {£-) <^ (£-) + {£+), and 
neglect all collisions involving two atoms in in either 
the final or initial state. The corresponding flux iV_|_ is 
thus of the form 

N+ = AiN^^ - A2N-N+. (1) 

The coefficients Ai and A2 are calculated using Boltz- 
mann equation |2^. In particular AiA^_ is exactly the 
evaporation rate in an evaporative cooling process [p5[. 
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FIG. 2: Integrated MCP signal during the RF probe pulse 
as a function of the delay between the truncation and probe 
pulses. The three curves correspond to 5 x 10^, 7 x 10^ and 
10 X 10^ atoms in the trap, varied by changing the power in 
the Zeeman slowing laser. The lifetime of the trap is 38 ± 4 s, 
and the temperature is 0.9 ±0.1 mK. 

If we make the further approximations that the atoms 
in £- and £+ have thermal distributions |p6| , neglect 
variations of the temperature during thermalization and 
assume that the colhsion cross section a is independent 
of velocity, Ai and A2 are analytic functions of the trap 
parameters, atomic mass to, a, rj and pL^Bo/kBT. This 
latter parameter appears because our trap cannot be 
approximated by an harmonic trap; we use the semi- 
linear form It is straightforward to take into ac- 
count the finite lifetime r of the atomic sample since 
N_{t) + N+{t) = iV_(0) exp (-t/r). The solution of the 
resulting differential equation is : 

l-g-exp[;^(l-e */^)] 

where T^h' = ^j^^ and iVt,, - (1 - g)A^-(O). 

The elastic collision rate is 701 = nav with n defined at 
the center of the trap and v = A^kBT/Trm. The quan- 
tities Vev, Ve and q are defined as in |2^, |2^; they are 
analytic functions of rj and ^3^0/^8^. The quantity q 
is the ratio of the number of atoms below the RF knife to 
the total for a thermal distribution (about 0.9 under our 
conditions), and Nth is the asymptotic value of for 
infinite trap lifetime. Numerical simulations of the en- 
ergy form of Boltzmann equation are in good agreement 
with our model for 77 > 10; for rj — 6 the quantity 7eiTth 
is 1.8 times larger meaning that our assumption about 
the distribution function fails for small 77 We take 
this factor into account in calculating 701. 

To fit the data of Fig. ^ with eq. (Q), we fix the lifetime 
T at its measured value and use Tth and A^th as adjustable 
parameters. The uncertainty in Tth is estimated by vary- 
ing the lifetime of the trap within its uncertainty range 
and looking at the resulting dispersion in Tth- The uncer- 



FIG. 3: Thermalization rate versus the area of the cor- 
responding TOF spectrum (proportional to the number of 
trapped atoms). The solid line shows a linear fit constrained 
to pass through the origin. 

tainty in the number of trapped atoms is estimated from 
the dispersion of the TOF area measurements before and 
after taking a curve as in Fig. ^. The exact value of q 
has little influence on the fit. 

We have made several tests to check the consistency of 
our results. First, we have checked that the fitted value of 
A'th corresponds to the expected fraction of atoms above 
the knife for our temperature. Second, Fig. |3| shows that 
^th^ is proportional to the number of trapped atoms, as 
it must be if the process of refilling of the upper energy 
class is due to two body collisions. We can exclude any 
effect independent of the number of atoms. The line pass- 
ing through the origin uses the slope as a fit parameter 
and has — 5 for 8 degrees of freedom. Third we have 
done an additional experiment that confirms the presence 
of elastic collisions: in a trap decay rate experiment, in 
the presence of the RF knife, the ion signal exhibits a 
clear non-exponential behaviour at short times. This ef- 
fect can be satisfactorily interpreted as elastic collisions 
bringing atoms above the RF knife and hence allows a 
measurement of the evaporation rate. This rate is con- 
sistent with the results obtained in our thermalization 
experiment. Fourth we have checked that heating can- 
not explain the repopulation of the upper energy classes. 
With the trap undisturbed, we can place an upper limit 
on the heating rate of 25 fiK in 60 s. This limit is two or- 
ders of magnitude too low to explain our data. Lastly, we 
have performed the thermalization experiment for differ- 
ent lifetimes of the magnetic trap (20, 40 and 60 s) and 
found consistent results. 

From our data in Fig. ^, we can deduce an accurate 
measurement of the thermalization time; the fit gives 
Tth — 3.0 ± 0.3 s for the densest sample. Using the mea- 
sured temperature and bias field, this value of Tth leads 
to 7oi = 6 ± 1 s^^; this result depends on the accuracy 
of our thermalization model. To find the rate constant 
a = jci/n, we must estimate the density. Since the data 
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show that our sample is close to thermal equilibrium, we 
can calculate the volume of the trap knowing the trap 
parameters. The absolute measurement of the number 
of atoms is performed by measuring the total power ab- 
sorbed from a saturating laser beam, similar to A 
TOF area of 75 mV.s corresponds to 10* atoms in the 
magnetic trap with an uncertainty of a factor of 2. This 
leads to a = 5 X 10~^ cm'^/s to within a factor 3 at 
T = 1 ± 0.1 mK. The ENS-Paris group has obtained a 
similar result with a different measurement p8[ . The uni- 
tarity limit at that temperature is a ~ 10~^ijin'^/s. This 
means that it is probably not valid to use a constant elas- 
tic cross section in our model and some deviation might 
appear in the quantity 7ciTth. We are currently investi- 
gating refinements to our thermalization model. 

The results shown here are very encouraging for 
evaporative cooling of He* in search of BEC. 
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